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ABSTRACT 

The peculiar velocities of galaxies are an inherently valuable cosmological probe, pro- 
viding an unbiased estimate of the distribution of matter on scales much larger than 
the depth of the survey. Much research interest has been motivated by the high dipole 
moment of our local peculiar velocity field, which suggests a large scale excess in 
the matter power spectrum, and can appear to be in some tension with the ACDM 
model. We use a composite catalogue of 4,537 peculiar velocity measurements with 
a characteristic depth of 33 /i _1 Mpc to estimate the matter power spectrum. We 
compare the constraints wi th this method, d irectl y studying the full peculiar velocity 
catalogue, to results from iMacaulav et alj (1201 II). studying minimum variance mo- 
ments of the velocity field, as calcu lated by I Watkins, Feldman fc Hudson! (|2009fl and 
iFeldman. Watkins fc Hudson] (|2010D . We find good agreement with the ACDM model 
on scales of fc > 0.01 h Mpc -1 . We find an excess of power on scales of fc < 0.01 h 
Mpc" 1 , although with a ler uncertainty which includes the ACDM model. We find 
that the uncertainty in the excess at these scales is larger than an alternative result 
studying only moments of the velocity field, which is due to the minimum variance 
weights used to calculate the moments. At small scales, we are able to clearly discrimi- 
nate between linear and nonlinear clustering in simulated peculiar velocity catalogues, 
and find some evidence (although less clear) for linear clustering in the real peculiar 
velocity data. 

Key words: cosmology: large scale structure of the universe - cosmology: observation 
- cosmology: theory - galaxies: kinematics and dynamics - galaxies: statistics 



1 INTRODUCTION & BACKGROUND 

The peculiar velocities of galaxies are a powerful cosmologi- 
cal probe, which directly traces the underlying dark mat- 
ter distribution (independent of galaxy bias) and is also 
sensitive to scales much larger than the size of the sur- 
vey. Recent interest in peculiar velocities has been driven 
by the high dipole moment of our local velocity field, 
which can appear to be in some tension with the ACDM 
model. The ACDM model indicates that we should ex- 
pect a peculiar velocity dipole of magnitude around 100 km 
s _1 , although many independent peculiar velocity surveys 
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show evidence for a bulk flow at low redshift of around 
400 km s -1 i n the direction I = 280 b = 10 degrees 
jHudson et ail 120041: I Watkins, Feldman fc Hudson! 2009: 



Feldman. Watkins fc Hudson 2010; Ma. Gordon fc Feldman! 



20111 ). However, iNusser fc David l|201ll ) find evidence for a 
flow more commensurate with ACDM of around 260 km s _1 
(in a similar direction as the other bulk flow results). 

In more tension with ACDM are significantly higher 
bulk flows at redshift ~0.25, which curiously appear to 
be in the same d i rectio n as the low redshift bulk flows. 
lAbate fc Feldman! (|2012l ) found evidence for an extremely 
high bulk flow of around 4000 km s _1 at redshift ~0.3, 
in a s imilar direction to other bulk flows. iKashlinskv et all 
j20ld ) find a bulk flow of around 1000 km s 1 extending to 
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z ~ 0.2 from kinematic Sunyaev-ZePdovich measurements. 
If substantiated, these bulk flows would indicate a non con- 
vergence of the peculiar velocity dipole, and present a seri- 
ous challenge to the assumptions of isotropy and homogene- 
ity. However, these bulk flow measurements are far from as 
robust as the directly measured peculiar velocity measure- 
ments we consider here; lAbate fc Feldmanl (J2012) note that 
their result may be entirely due to systematic effects. In 
this work, we will only consider the direct peculiar velocity 
measurements at low redshift. 

In addition to the high dipole moment, there also ap- 
pears to be a low shear of the velocity fie ld (|Jaffe fc Kaiserl 
1 19951 ; iFeldman. Watkins fc Hudsonl |2010| ). which indicates 
that the density contrast responsible for the velocity dipole 
is on extremely large scales. At the depths of the peculiar 
velocity surveys (up to 100 /i _1 Mpc), this suggests an excess 
density contrast on scales ~ 1 h~ 1 Gpc. Th e volume of space 
probed by gala xy redshift surveys (e.g. ICole et al 1 120051 ; 
iReid et "ali i 2010 and Blake et al. 201~cj) is currently too small 
to robustly constrain clustering on these scales, although 
Prii omas, Abdalla fc Lahavl i|201ll ) found evidence for excess 
large scale power in the MegaZ photometric redshift survey. 
Measureme nts of the CMB can probe anisotropics on these 
large scales |Hlozek et al.ll2012l ). although to compare to re- 
sults at low redshift, the growth of these anisotropics must 
be assumed, which depends on the cosmological model. 

There are many possible explanations for the high 
dipole moment, ranging from systematic effects to more 
exotic explan ations invoking exten ded cosmology or mod- 
ified gravity. iHudson et al. l|2004 ) considered systematic 
effects in the SMAC (Streaming Motions of Abell Clus- 
ters) peculiar velocity survey, such as a dipole variation 
in the velocity dispersion, galactic extinction, and calibra- 
tion across different observations. They found that sys- 
tematic effects could account for at most half of the high 
dipole moment of the SMAC survey, leaving a dipole 
moment which is still at least three times higher than 
the ACDM expectation. There are many theoretical pos- 
sibilities to extend the ACDM model to produce a large 
scale excess of power and a p eculiar velocity dipole, such 
as m o difications to grav it y. (Avaita. Weber fc Wetterichl 
120091 : 



.Khourv fc Wvmanl 120091) dark energy cluster- 
ing l|Potter fc Chongchitnanl l201l]'). vorticitv. dPalld I201C 1 



or 'tilted' universes ([Mcrsini-Houghton & Holman 200E 
iKashlinskv et al.1 12010| ; iMa. Gordon fc Feldmanl |201lf ). In 
this work we do not address any theoretical explanation 
in particular, but by using the peculiar velocity data to 
constrain the power spectrum in model independent band- 
powers we aim to provide results which will be useful to 
constrain a range of explanations for the high dipole. 

Jaffe fc Kaiser] (Il995h IKolatt et ail (1 19961) 



IZaroubi et al l ll 19971). IKolatt fc Dekell (|l997l ). IZaroubi et al.1 



( 200ll ) and ISilberman et al.l (|2001 ) have used peculiar 



velocity surveys to infer the matter power spectrum directly 
at z = 0. More recently in iMacaulav et all l|201ll ). we in- 
ferred the underlying power s pectrum from moments of the 
peculi ar v elocity field from Watkins, Fel dman fc Hudson! 
(|200gh and lFeldman. Watkins fc Hudsonl (|2010l ). specifically 
in the context of understanding the high bulk flow. One of 
our main findings was that the excess of power indicated 
by the anomalously high dipole moment was dramatically 
reduced when the shear and octupole moments were 



also included. This leads us to ask if the correspondence 
with the ACDM model would be improved if we were to 
hypothetically include higher still moments of the velocity 
field. That is the motiva tion for this work, a lthough we 
take the approach used by I Jaffe fc Kaiserl (|l995l ) to analyse 
peculiar velocity catalogues directly, without compressing 

the data into moments. 

Recently, lAbate fc Erdogdul (|2009h applied a similar 
formalis m to constrain the modif i ed gr avity parameter 7. 
Similarly IMa. Gordon fc Feldmanl (|201ll ) used a similar for- 
malism to fit for parameters of the dipole moment and ve- 
locity dispersion parameter. In both fiducial ACDM 
power spectrum was assumed, and additional parameters 
of interest were allowed to vary. The key difference between 
those papers and this work is that here we treat the underly- 
ing power spectrum as a set of free parameters, as opposed to 
fitting extra parameters for additional effects beyond a fixed 
power spectrum. In this way, we obtain new constraints on 
the power spectrum which are independent of the fiducial 
cosmology. 



2 METHOD 

In this paper we consider two distinct methods to relate 
peculiar velocity measurements to large scale structure: A 
maximum likelihood based approach to analyse a full pe- 
culiar velocity catalogue (the 'catalogue' method), and an 
alternative approach studying minimum-variance moments 
of the velocity field (the 'moments' method). 

As well as an apparent radial velocity due to the Hub- 
ble flow, galaxies also have a peculiar velocity towards local 
over-densities of matter. This peculiar velocity field v(r) can 
be related to the matter density contrast S by 



4tt 



r' 



(1) 



where f g is the growth rate of th e density con trast, 
9m<5/<91na, and a is the scale factor |Peeblesl fr993). The 
density contrast is defined in terms of the ratio of the den- 
sity at r, p, to the average density p, so that 8 = p/p. 

We can measure the peculiar velocity via the effect it 
has on the redshift of the galaxy. The redshift z of a galaxy 
is given by a contribution from the Hubble flow, Hor, and 
the line of sight component of the peculiar velocity S: 



cz = H r + S 



(2) 



where, for a galaxy labelled m, at position r m , we have the 
line of sight peculiar velocity, S m given by 



(3) 



where r m is a unit vector in the direction of galaxy m. 
To measure S we thus need to combine the redshift of 
the galaxy with an independent measure of the distance r. 
In principle, we can calculate the peculiar velocity of any 
galaxy for which we have the luminosity distance, measured 
using distanc e indicators such as supernovae (SN), Tully- 
Fishe r (TF) (iTullv fc Fisherl Il977l ) or Fundamental-Plane 
(FP) (|Dressler et al.lll987t ) measurements. Individual uncer- 
tainties on luminosity distances are typically rather large 
(5% for SN, and around 10 to 20% for TF and FP), which 
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propagates to a very large uncertainty in the peculiar veloc- 
ity. 



2.1 Maximum Likelihood Catalogue Method 

We now consider how to relate a catalogue of peculiar veloc- 
ity measurements S m t o large scale structure , following the 
method as presented in ljaffe fc Kaiserl ()l995l ). The method 
is based on a likelihood framework, where the likelihood C 
is given by 



2n N / 2 \R„ 



[1/2 eX P ^SrnRmnSr 



(4) 



where Rmn is the covariance matrix for the peculiar velocity 
measurements, for a catalogue of TV galaxies. The covariance 
matrix is split into two components, a 'velocity' and an 'er- 
ror' term: 



Jin 



(•») 



(v) (e) 

Rmn models the coherent large scale structure, and R mn is a 
noise component to account for nonlinear velocity dispersion 
<t*, and uncertainty in each peculiar velocity measurement 
<r m . If we assume that the measurement errors are uncor- 
related, Rm n simply contributes these uncertainties to the 
diagonal of the covariance matrix, and is given by 



Rmn = (d"m + &*)5„ 



(6) 



where S mn is a Kronecker delta. This formalism does not 
include correlations between measurement uncertainties. In 
principle, there may be small correlations between measure- 
ment uncertainties due to, for example, galactic extinction, 
or the assumed mass-to-light ratio within catalogues of the 
same type of galaxy. However, we assume that the cross- 
correlations in the data covariance matrix are dominated by 
the geometry of the survey, which we model with the 've- 
locity' component of the covariance matrix. This term of 
the covariance matrix contains the power spectrum we are 
fitting for, and is given by 



Ji-mn (^0 



ink dk 
(2tt) 3 



Pv (&),/*mn (&) 



where P v (k) is the velocity power spectrum, which is related 
to the matter power spectrum by 



(8) 



Here a is the scale factor. The window function f mn (k) is 
given by 



fmnik^j — ^m,i^"n,j 



d 2 k 
4% 



ik\z- (r m — r,, ) 



(9) 



and can be calculated analytically in terms of trigono- 
metric functions. A derivatio n of f mn is presented in 
iMa. Gordon fc Feldmanl (|201lh . The general approach we 
take is to map out the likelihood in terms of parameters of 
the power spectrum. We next consider an analogous method 
to analyse moments of the velocity field, before considering 
parametrisation of the power spectrum any further, since 
the parametrisation is common to both methods. 



2.2 Minimum Variance Moments Method 

In this paper we directly compare results between the maxi- 
mum likelihood catalogue method to an alt ernative method 
study ing moments of the velocity field, from lMacaulav et al.l 
(|20 1 IT ) . We can consider the velocity field as a Taylor expan- 
sion, given by: 



Vi(r) =Ui + UijTj + U l j k r ] r k + 



(10) 



where Ui is the dipole moment of the velocity field (often 
called the 'bulk flow'), and provides most information about 
the largest scale fluctuations. Uij is the shear of the velocity 
field, sensitive to intermediate scales. Uijk is the octupole 
moment, and is sensitive to scales less than the size of the 
survey. The relative sensitivity of each moment is shown in 
Figure [T] Since we can only measure the line of sight compo- 
nent of each peculiar velocity, each individual measurement 
must be weighted acc ording to the component of t h e mo - 
ment it is sensitive to. Watkins. Feldman fc Hudson! (|2009h 



and iFeldman. Watkins fc Hudson! (|2oiof T" developed new 
'minimum variance' weights to estimate the velocity of the 
volume traced by the galaxies in the survey. These weights 
are designed to minimise the effects of small scale motions, 
to provide a better estimate of the large scale velocities. We 
do not reproduce the derivations of the weights he r e; the 
method is presented in lWatkins. Feldman fc Hudson! (I2009T) 
for th e dipole and extended in lFeldman. Watkins fc Hudson! 
(2010) for the shear and octupole. 

With this method, the data consists of the three com- 
ponents of the dipole vector, six independent components 
of the shear, and ten independent components of the oc- 
tupole. We take a similar approach as with the catalogue 
method, splitting the covariance matrix into 'velocity' and 
'error' terms. The 'velocity' term is now given by 



R$ = I dkP{k)W 2 m {k) 



(11) 



where p and q index the 19 independent moments. The win- 
dow function Wp 9 (fc) is sensitive to different scales for the 
dipole, shear and octupole, as plotted in Figure [1] We con- 
struct a likelihood in exactly the same manner as the cata- 
logue method, in terms of parameters of a power spectrum, 
which we now consider. 



2.3 Power Spectrum Parametrisation 

Ultimately, we wish to constrain parameters of the under- 
lying matter power spectrum. There are many choices for 
ways to parametrise the power spectrum. A popular choice 
in many earlier works was the power spectrum shape pa- 
rameter F, the matter density Q. m and the power spec- 



trum normalization a t 8 h 1 Mpc. cr» |jaffe fc Kaiserl Il995l ; 
iKolatt fc DekelHl995l ). ISilberman et all (|200lh opted for a 



band-power parametrisation, which is the approach we take 
here, in terms of band-powers, P a , so the power spectrum is 
given by 



P(k) = 



Pa ka < k < k a +i 
otherwise. 



(12) 



The velocity covariance matrix is then 
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Figure 1. Comparing sensitivity to different scales for the mo- 
ments and catalogue methods. For the moments panel, we plot 
the average of the diagonal of the covariance matrix for the dipole, 
shear and octupole terms, for COMPOSITE. Since the moments 
are designed to be orthogonal, the diagonal of the covariance ma- 
trix gives an indication of the scales at which the moments are 
most sensitive. For the catalogue panel, the diagonal of the covari- 
ance matrix is 1/3, so we plot the sum of the entire covariance 
matrix, normalised by the square of the number of galaxies in 
each catalogue, which illustrates the scales at which the window 
function is sensitive. The sub-catalogues which span the greatest 
distance, DEEP and SMAC, are most sensitive to largest scales, 
while the shallowest catalogues, SFH — h 9 and ENEAR, are most 
sensitive to smaller scales. 



where K, is given by 



(13) 



(14) 



This parametrisation is directly sensitive to the combination 
of fgP a - To directly constrain P{k), we must either assume a 
fiducial growth rate, or marginalise over other measurements 
of the growth rate. For generality, we choose to treat the 
combin ation of f a P a as a para meter. 

In iMacaulav et all (|201lD . we used flat band-powers. 
This parametrisation was sufficient to demonstrate the large 
effect of including just the dipole, or higher moments of the 
velocity field. However, in the widest, single band parametri- 
sation, the flat bands can introduce an apparent artificial 
shift towards an excess of power, shown in the upper panel 
of Figure[2] We find that when we factor in a fiducial ACDM 
power spectrum into the window function, P^(k), and allow 
a constant amplitude of this power spectrum, A a , to vary, we 
obtain an amplitude which is consistent with ACDM. That 
is, for the same average amplitude of P(k) the most likely 
value can depend on the slope of the band within the band 
range. We find in practice that this shift is most significant 
for the single band parametrisation, and when we allow sev- 
eral bands to vary across the range the slope of each band is 
less significant. Nevertheless, for consistency and complete- 
ness we use ACDM shaped bands for all parametrisations, 
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Figure 2. The effect of band shape. In the upper Figure, when 
using the constant amplitude of a single flat band-power to 
parametrise the power spectrum, we find a systematic excess com- 
pared to the average ACDM power spectrum in the same k range. 
However, when we use the amplitude of the ACDM power spec- 
trum, A a as a parameter, we find no systematic deviation from 
the A a = 1 expectation. The results are for five different realisa- 
tions of a ACDM simulation, each consisting of dipole, shear and 
octupole moments. 
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Figure 3. One standard deviation contours for one power spec- 
trum band and the velocity dispersion <r„ , for moments generated 
from a simulated velocity field. The contours are fairly insensi- 
tive to the choice of cr*. Since compressing the velocity field into 
moments is insensitive to small scale effects, these contours are 
consistent with c* = 0. The Id likelihood distributions are for a 
fixed (T* of 200 kms" 1 



varying the average amplitude within each band. The kernel 
in this case is 



ICa 



dkf mn (k)P M (k) 



so our velocity covariance matrix is now 
(Hf 



R 
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where our parameter is now A a , the amplitude of the fiducial 
power spectrum in each band. We thus expect A a — 1 in all 
bands for a ACDM velocity field. To more easily compare 
A a to expectations, we multiply the results by the mean of 
the power spectrum in each band, {P a ), given by 



550 



(17) 



Although this parametrisation specifies the shape of the 
power spectrum within each band, we see in Section [4] in 
simulated catalogues that the method is fairly robust to the 
choice of fiducial power spectrum. We find that the cata- 
logue method is sensitive to the velocity dispersion param- 
et er, and thus treat it as a free parameter in the manner 
of iMa. Gordon fc Feldmanl (|201ll ). To directly compare be- 
tween the catalog ue and moments metho d, we repeat much 
of the analysis of iMacaulav et all (|201lT > with a, as a free 
parameter, although we find with results from velocity mo- 
ments that marginalising over the velocity dispersion param- 
eter has only a very small effect. The same set of moments 
from Figure [2] are analysed in Figure allowing both the 
velocity dispersion and amplitude to vary. In this paper we 
now consider the combination of fgP a as a parameter for 
both met hods, as opposed to m arginalising over the growth 
rate as in IMacaulav et all (|201ll h We assume a fiducial value 
for the velocity dispersion of 350 km s _1 . 

We use a Metropolis-Hastings Markov-Chain Monte- 
Carlo method to map out the likelihood expression in 
terms of the band-powers and velocity dispersion. The er- 
ror bars plotted on the power spectrum plots are calcu- 
lated at the one standard deviation confidence level from the 
marginalised likelihood distribution for each band-power. 
Where the likelihood peaks at zero this is plotted without a 
marker, and the error bar represents the extent of the upper 
bound. 



3 PECULIAR VELOCITY DATA 

We primarily study the COMPOSITE peculiar velocity 
catalogue, which consists of 4537 individual peculiar ve- 
locity measurements, and a characteristic depth of 3 3 
/i _1 Mpc, compiled by iFeldman. Watkins fc Hudson! l)2010h . 
The COMPOSITE catalogue is composed of several sub- 
catalogues, which we also consider individually. The largest 
sub-catalogue is the SFIH — h (Spiral Field I-band) sam- 
ple, which consists of 3456 TF measurements. We anal- 
yse the field galaxies and groups in the SFI++ sam- 
ple separately, as SF I++_£ (2720 galaxie s ) and SFI+- |-o 
(736 m easurements) l|Masters et alj 120061 : ISpringob et al.l 
2009a |3). We also study the combined DEEP ca talogue, 
compiled by IWatkins. Feldman fc Hudson! l|2009l ). which 
consists of 294 of the deepest peculiar velocity measurements 
and the ENEAR (Ea rly-type NEARby galaxies) catalogue 
ijda Costa et al.ll2000h . The smallest c atalogue we conside r 
individually is the SMAC catalogue l|Hudson et al.ll2004h . 
which is a sub-set of the DEEP catalogue. 

We test our procedures on peculiar velocity catalogues 
generated from ACDM simulations. We analyse the set of 20 
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Figure 5. Constraints on the velocity dispersion parameter a, 
and one band-power, in the range k = 0.002 to 0.196 h _1 Mpc, 
for simulated peculiar velocity surveys from the Las Damas and 
Mark III simulations. The markers represent the peak of the like- 
lihood for each catalogue. The contours are 1 standard deviation 
uncertainty about the peak likelihood, which we have only plot- 
ted for three catalogues for each simulation set, to prevent the 
plot from becoming overcrowded. Similarly, we only plot results 
for 10 of the Mark III catalogues. These uncertainty ranges are 
typical for the points plotted here. 



simulated Mark III c atalogue^] drawn from realisations of 
the Virgo simulation ()Kolatt et al.l!l99r3 ). consisting of 1300 
entries each (some of the catalogues consisted of slightly 
more than 1300 entries - these were trimmed to 1300 for 
consistency). We also study six simulated catalogues from 
the Las Damas simulation (McBride et al., in prep.), de- 
signed to resemble the SFI++ / catalogue with 2720 entries. 
Histograms of the real and simulated surveys are shown in 
Figure [4] 

We study the simulated data sets with both the cata- 
logue and moments method. We analyse the real data with 
the catalogue method for the combined COMPOSITE set, 
and the five sub-catalogues. We find that the sub-catalogues 
are too small to reliably study the moments alone, so we 
present results here with the moments method only for the 
COMPOSITE catalogue. 



4 RESULTS 

4.1 Simulated Data 

We begin by testing our method on peculiar velocity cat- 
alogues generated from ACDM simulations. We start the 
analysis with the simplest parametrisation: one band-power. 
We choose a window function range of k — 0.002 to 0.196 
fe -1 M pc, to match the k range used in IMacaulav et alj 
l|201lf) . We also include the velocity dispersion as a free pa- 
rameter. Results are shown in Figure [5] 

We find that for this band-power parametrisation, the 
velocity dispersion is uncorrelated with the power spectrum 
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Figure 4. Histograms of the real surveys, and the Las Damas and Mark III simulated surveys, illustrating the distribution of line of 
sight peculiar velocities, 5, uncertainty in S, erg, distance r, and coordinates I and b. The COMPOSITE survey is plotted as a stacked 
histogram and is comprised of the SFI++/ (black), SFI++ 9 (green), ENEAR (purple), DEEP (red), and SMAC (brown) surveys. 



amplitude. The velocity dispersion is slightly higher than 
the ACDM expectation for the Las Damas mocks, although 
it is lower in the Mark III mocks. In both sets of mocks, 
we note a small systematic shift from the ACDM expecta- 
tion. As this shift is in opposite directions for the Mark III 
and Las Damas mocks, it seems reasonable to conclude that 
this may be due to the many particularities of generating 
a mock catalogue from an underlying power spectrum, as 
opposed to a systematic shift introduced by our method. 
In other words, the Mark III and Las Damas mocks taken 
together surround the ACDM expectation well. Following 
iMacaulav et al.l (|20Tlf ). we next consider a parametrisation 
with three band-powers, spanning the same k range with 
bands evenly spaced in log k. The results are shown in Fig- 
ure [5] 

We will refer to the bands from the largest scales to the 
smallest scales as bands 1 to 3. We find that the uncertainty 
in band 1 is larger in the simulated data with the catalogue 
method than the moments method (not plotted). The uncer- 
tainty in band 1 is particularly important for understanding 
the high dipole moment, and will be discussed further in 
Section [5] Both sets of catalogues are much more sensitive 
in band 3 than with the moments method; we thus decide 
to extend that band-power parametrisation by adding two 
more bands (4 and 5) at smaller scales. Results are shown 
in Figure [7] 

We find that band 5 is anti-correlated with the velocity 
dispersion parameter. This is not surprising, given the ~ 5 
/i _1 Mpc scales spanned by this band. We are encouraged 
to see a difference at over la from the nonlinear Halofit 
l|Smith et al.ll2003T ) corrections to the galaxy power spectra, 
illustrating the ability of peculiar velocity measurements to 
directly probe the underlying linear matter distribution. As 
a test, we repeat the results in Figure [7] assuming a fiducial 
Halofit corrected power spectrum to factor into the window 
function. The results are shown in Figure [5] Encouragingly, 
we see in the smallest two bands that the data still favour 
the linear power spectrum. 




Figure 6. Results from the Mark III and Las Damas mocks in 
three band-powers. The velocity dispersion was varied as a free 
parameter, and marginalized over in these results. The grey curve 
is the ACDM power spectrum. The horizontal black lines span 
the k ranges of the band-powers, and indicate the average value 
of the power spectrum across this range. The results should be 
compared directly to these. The markers are the marginalized 
results for each band-power, for each of the Mark III and Las 
Damas surveys. The k location within each band is arbitrary; 
each point should be fully considered equally at the centre of the 
band, spanning the full width. 



4.2 Real Data 

We now consider results from real peculiar velocity surveys. 
We start with the one band-power and velocity dispersion 
parametrisation, shown in Figure [5] 

Using the ACDM shape band, we find a slightly 
lower average amplitu de than the flat band result from 
IMacaulav et alj l|201ll ). Although we note a small excess, 
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Figure 7. As Figure [6] but with two additional bands added 
at smaller scales. The smallest scale band is at scales for which 
galaxy power spectra must be corrected for (e.g., with Halofit). 
With the peculiar velocity method, we are sensitive to the total 
matter distribution, and can thus directly probe the linear power 
spectrum, without accounting for halo corrections. 
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Figure 9. 1 standard deviation contours from the real data, 
parametrised in terms of the velocity dispersion ct* and one band- 
power, in the range k = 0.002 to 0.196 /i _1 Mpc. There is a slight 
excess of power in all the catalogues, although within the la level 
for most. 
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Figure 8. As Figure [7] but with a fiducial power spectrum cor- 
rected with Halofit. The black bars now represent the mean value 
of the Halofit fiducial power spectrum which was factored into 
the window function. The data still favor the linear power spec- 
trum, indicating that the method is robust to the choice of fiducial 
power spectrum. 



the ACDM value is now well within the la uncertainty 
range of the moments result. We find that the results for 
the velocity dispersion agree w ell with similar results from 
iMa. Gordon fc Feldmanl (|201lf ) . We find a smaller uncer- 
tainty in the band-power with the catalogue method than 
the moments method, and a much improved constraint on 
the velocity dispersion parameter. We next consider a three 
band parametrisation, shown in Figure [UJl 

With the moments method, band 3 was merely an upper 
limit - it is well constrained with the catalogue method. The 
low shear of the velocity field caused band 2 to be lower 
than the ACDM expectation. When we now analyse the full 




Figure 10. As Figure|6] for the real peculiar velocity catalogues. 
As before, the k position of each point within the bin is arbitrary, 
and each point should be fully considered at the centre of each 
band. The points have been ordered in each band from left to right 
by the number of galaxies in each catalogue. We also plot results 
from the moments analysis of the COMPOSITE catalogue. 



catalogue, the band agrees extremely well with the ACDM 
expectation. We still observe an excess of power in band 
1, although the lower bound is not constrained here. These 
results are discussed further in Section [5] We next extend 
the parametrisation with two smaller bands, as before with 
the simulated catalogues. The results are plotted in Figure 
□3] and presented for COMPOSITE in Table [U We find 
good agreement with ACDM in bands 3 and 4. We find that 
bands 1 and 2 are anti-correlated. 
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Figure 11. As Figure [7] now for real data. In band 1 we see 
a slight excess of power, although the uncertainty includes the 
ACDM model at the lcr level. Band 5 favours the linear clustering 
power spectrum over the nonlinear model. 
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/fP(fe) [(^Mpc) 3 ] 
ACDM COMPOSITE 



4.2 x 10 3 
6.2 x 10 3 
2.2 x 10 3 
2.8 x 10 2 
1.6 x 10 1 



2.5" 1 



l-9li 9 x 10 4 
x 10 3 



3.0 4 



0.5 



+5.7 
-0.5 



X 10 1 



Table 1. Results for the COMPOSITE catalogue five-band 
parametrisation, as plotted in Figurc fTTI Although the high dipole 
velocity leads to an excess of power by over a factor of 4 in band 
1, the uncertainty includes the ACDM model within the la un- 
certainty level. 



4-2.1 Comparison to \Silberman et al\ 1(200 A) 

We also test our metho dology by reproducing the results of 
ISilberman et al.l (|200lf ). to which our method is similar. We 
change our bands as follows, to match the bands used in their 
work: For k ^ 0.02 ftMpc -1 we use a ACDM band, which we 
do not vary. In the range 0.02 < k 0.07 and 0.07 < k sC 0.2 
/iMpc" 1 we use two independent flat bands of constant am- 
plitude. For k > 0.02 /iMpc^ 1 we use a power law band, 
A a k n The spectral slope n is set to -0.95 for the simulated 
Mark III catalogues, and to -1.4 for the real SFlH — h/ cata- 
logue, and the free parameter A a is varied. The results are 
shown in Figure 1121 We rep roduce the slight decre ment at 
k ~ 0.1 ftMpc -1 as noted bv lSilberman et al.l (|200lT ). which 
we also observe in the SFI++ / catalogue in band 3 of our 5 
band parametrisation in Figure QT] 



4.3 Fitting a Velocity Dipole 

In the three band model, we find that the uncertainty in 
band 1 is larger with the catalogue method than with the 
moments method. We see these results in the simulated cat- 
alogues, but it is particularly important to understand the 
result in the COMPOSITE catalogue since it is closely re- 
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Figure 12. Analy zing a simulated Mark H I catalogues with the 
parametrisation of :S ilberman et all |200l]) . Th e results we find 
here a gree well with the results in Figure 7 of of Silbcrman ct al. 
J200lh . 



lated to the high velocity dipole. As can be seen in Figure [T] 
the catalogue method window function is most similar to the 
window function of the shear moment, and (relatively) not 
as sensitive to the scales probed by band 1. To try and under- 
stand these results, we use our maximum likelihood method 
to estimate t he dipole, to compare to the minim um vari- 
ance dipole of iFeldman. Watkins fc Hudson! ()2010l ). Instead 
of varying the power spectrum, we fix the power spectrum 
at the fiducial ACDM value and minimise parameters o f the 
dipole, as analysed bv lMa. Gordon fc Feldmanl (|201tJ ). We 
find that it is not possible to simultaneously constrain the 
thr ee band-power spectrum and a velocity dipole. Follow- 
ing [MaZGordon^^idma^ l|201ll 'l. we model the dipole U 
by subtracting the line of sight component of the dipole to 
each galaxy in the catalogue from it's line of sight peculiar 
velocity, S n , to obtain a 'tilted' velocity, p n 



Pn — S n 



u 



(18) 



We parametrise the dipole bulk flow as U = {U r , Ui, Ub}, 
where U r is the magnitude of the bulk flow (in km s" 1 ), and 
Ui and Ub are the direction of the flow in degrees (galactic 
coordinates). We also include the velocity dispersion <t» as a 
free parameter, which we marginalise over to find the dipole 
bulk flow. We find a dipole of mag nitude U r = 380±l 9 i2 
km s _1 in the direction U t = 295±^ and U b = 14 ± 18 
degrees. The direction of the dipole is shown in Figure 1131 
This is slightly slower and w ith a larger uncertainty than 
the di pole moment found by IFeldman. Watkins fc Hudson! 
(|2010h of U r = 416 ± 78 km s _1 in the direction Ui = 282 ± 
11 and Ub = 6 ± 6 degrees, and is related to the larger 
uncertainty we find in band 1 with the catalogue method 
than the moments method. A comparison of measurements 
of the dipole is shown in Table [2] 

Results from many different surveys and analysis meth- 
ods appear to agree on the direction of the velocity dipole, 
at around I — 280 and b = 10 degrees, within an uncertainty 
radius of about 15 degrees. However, there is less consensus 
as to the magnitude of the dipole. The magnitude of the 
dipole (and the corresponding uncertainty) can depend on 
the depth of the survey, the treatment of outliers, and the 
weighting of galaxies. Even so, we can very conservatively 
state that the magnitude of the velocity dipole out to 100 
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U r [km s" 



-1 



Ui [degrees] Ub [degrees] 



Watkins. Feldman 



Nusscr 




Davis (2011) (SFI++, 40 fc^Mpc ) 
(SFI++, 100 fc-iMpc ) 
This work 



407 ± 81 
416 ± 78 
340 ± 130 
333 ± 38 
257 ± 44 
380 4 



-,+99 
-132 



287 ± 9 
282 ± 11 
285.1+f 9 ^ 
276 ± 3 
279 ± 6 
295± 18 



8 ± 6 

6 ± 6 

a 1 +is.5 

a -17.8 

14 ± 3 
10 ± 6 
14 ±18 



Table 2. 



Comp, 



arison of dipole 



measurements within 100 h x Mpc. This work, Watkins. Feld man &; H udson! ll2009t), 
IFeldman, Watkins &: Hudson! ll20iri ) and lMa, Gordon k, Feldman! J201lf) all study the COMPOSITE catalogue, whereas iNusser fc Davisl 
j201lh study the SFIH — h catalogue, which comprises the majority of COMPOSITE. The direction of the flow from different works agree 
well, although there is considerable variation in the magnitude of the flow. The magnitude of the flow can depend strongly on the depth 
of the survey, how galaxies at different depths are weighted, and the sensitivity of the method the fiducial power spectrum. 




COMPOSITE 

Ma, Gordon & Feldman 20 1 1 
Nusser& Davis 2011 
Feldman, Watkins & Hudson 2010 



Figure 13. The direction of the best-fitting velocity dipole 
in the catalogues, (in a Mollweide projection of the sky 
in galactic coordinates). The contours represent 1 standard 
deviation levels on the dipole direction angles Ui and Ui, 
found in the peculiar vel ocity catalogues. The directio n agree s 
with other results from IFeldman, Watkins fc Hudson! 201fj|) , 
iMa. Gordon fc Feldmar] l|201lh and lNusser fc Davisl (20111 )! 



h Mpc is a factor of several times higher than the ACDM 
expectation. 



5 DISCUSSION & CONCLUSIONS 

In terms of understanding the high dipole moment, the main 
difference between the two methods we have considered here 
is that the catalogue method finds an excess of power in band 
1 with an uncertainty which includes the ACDM model at 
the lcr level, while the moments method finds an excess of 
power with an uncertainty which excludes the ACDM model 
at over the la level. This is the issue we now consider. The 
power spectrum at these scales is closely linked to the mag- 
nitude of the dipole moment, where we similarly find a larger 
uncertainty with the catalogue method than the minimum 
variance moments result. The key to the different results 
are the minimum variance weights, which are designed to 
be sensitive to the largest scales, at the expense of small 
scale information. Indeed, of the sub-catalogues we analyse 
with the catalogue method, band 1 is best constrained by 
the DEEP and SMAC samples, which are the deepest of 
the sub-catalogues considered here. While both DEEP and 
SMAC are included within COMPOSITE, we find that the 
additional, shallower galaxies are effectively acting as noise 



as far as constraints on band 1 are concerned. Essentially, 
the minimum variance weights achieve this effect to a maxi- 
mal extent: preferentially weighting the deeper galaxies with 
the cleanest measurement of the velocity moments, at the 
expense of small scale information. 

On smaller scales, we also find a difference between the 
two methods in band 2. We find that the catalogue method 
agrees extremely well with the ACDM model, while the mo- 
ments method underestimates the power. This is due to the 
low shear of the velocity field, which provides most sensitiv- 
ity at the scales of band 2. When we analyse the full velocity 
catalogue, the small scale motions not modelled by the mo- 
ments appear to combine to provide a constraint which again 
agrees extremely well with ACDM. In much the same way 
as the one band parametrisation in iMacaulav" et all (|2011l ) 
illustrated the difference of including only the dipole, or ad- 
ditionally the shear and octupole, the result in band 2 high- 
lights the effect of considering only moments of the velocity 
field, or the full information available to us. 

In the parametrisation we have chosen, we have im- 
posed that P(k) = outside the range of the band-powers, 
which could, in principle, affect the results in the bands at 
the edge of the parametrisation. However, for the five band 
parametrisation, the range of the bands is sufficiently wide 
that the entire sensitivity range of the catalogues is cov- 
ered, and the bands are insensitive to power beyond these 
scales. We can also directly assess this effect by comparing 
the three and five band models we have presented here. In 
the three band model, we have that P(k) = at the scales 
where P(k) 7^ in the five band model (when we include 
bands 4 and 5). However, we find that the results in band 
3 are similar for both models (as can be seen in figures 
[7] for the simulated catalogues, and Figures [T7J1 fc ITT1 for the 
real catalogues), whether the neighbouring band is zero (the 
three band model) or nonzero (the five band model). 

We can also consider the issue of whether our method 
is robust to systematic errors in the peculiar velocity data. 
For example, the Tully-Fisher and Fundamental Plane mea- 
surements we consider both depend on the assumed mass- 
to-light ratio of the observed galaxies. To quantify the effect 
of systematic errors in the distance indicators, we have re- 
analysed a simulated Mark III catalogue and the genuine 
SGI++ 3 catalogue, where we have artificially introduced a 
systematic offset in the distance indicator and then repeated 
the power spectrum analysis. This mimics a systematic er- 
ror which may be introduced by incorrectly determining the 
intrinsic luminosity of the galaxies. We have analysed each 
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catalogue with eight different artificial offsets in the distance 
indicators in the range from an underestimate of the dis- 
tance by a factor of 0.6 (corresponding to an overestimate 
of the peculiar velocity) to an overestimate of the distance 
uncertainty by a factor of 1.4 (corresponding to an underes- 
timate of the peculiar velocity. As expected, we find a con- 
sistent overestimate of the band-powers when the peculiar 
velocity is overestimated, and a consistent underestimate 
in the band-powers when the peculiar velocity is underes- 
timated. We find that the ultimate systematic error in the 
bands scales linearly as approximately twice the error of the 
artificial distance offset. 

However, a robust indication that such systematic er- 
rors in the distance measurements are not problematic is 
that we observe similar results in each of the sub-catalogues 
of COMPOSITE that we study. The DEEP and SFI++ 
catalogues consist of different distance indicators, and we 
observe similar results both for the power spectrum fitting 
and dipole fitting; there is no individual catalogue which is 
anomalous. 

Perhaps more problematic may be systematic errors 
which can introduce a dipole variation in the observed lu- 
minosities, such as galactic extinction. Such a systematic 
effect would affect all the catalogues in the same manner, 
and could also mimic a bulk flow. From analysing moments 
of the velocity field, we know that a velocity dipole has most 
effect on the large scales of bands 1 and 2; the results we find 
in bands 3, 4 and 5 should be fairly robust to any systematic 
errors which introduce a dipole variation. As for bands 1 and 
2, although some of the velocity dipole could, in principle, 
be due to galactic extinction, it would then be difficult to 
reconcile with the low shear of the velocity field, and even 
harder to mimic with systematic effects. In other words, the 
high dipole and low shear, when considered together, are in- 
dicative of the large scale over densities we see in band 1, 
and are difficult to fully mimic with systematic effects. 

To conclude, we find that inferring the underlying power 
spectrum from peculiar velocity cata logues continues the 
general trend of lMacaulav et all (|201ll ) : that including more 
detail in the velocity field improves the agreement with the 
ACDM model. Specifically, we observe good agreement with 
ACDM on scales of k > 0.01 /iMpc -1 , although the agree- 
ment with ACDM in band 1 from the catalogue method 
is only due to the larger uncertainty than the moments 
method. While the high dipole moment alone may appear 
anomalous, when we consider the full peculiar velocity mea- 
surements, we find a power spectrum which agrees well with 
the ACDM model. 
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